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Abstract 

We examine Bose-Einstein condensation as a form of symmetry breaking in 
the specific model of the Einstein static universe. We show that symmetry 
breaking never occurs in the sense that the chemical potential fi never reaches 
its critical value. This leads us to some statements about spaces of finite 
volume in general. In an appendix we clarify the relationship between the 
standard statistical mechanical approaches and the field theory method using 
zeta functions. 
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I. INTRODUCTION 



Bose-Einstein condensation (BEC) for non-relativistic spin-0 particles is standard text- 
book material In the infinite volume limit, there is a critical temperature at which a 
phase transition occurs. For a real system, such as liquid helium, the effects of interactions 
may be important. (See [Q for a recent review.) The study of Bose-Einstein condensation 
for relativistic bosons is more recent. In particular, Refs. [0-0] applied the methods of rel- 
ativistic quantum field theory at finite temperature and density to study BEC. The phase 
transition, which occurs at high temperatures, can be interpreted as spontaneous symmetry 
breaking. Subsequent work extended the analysis to self-interactions in scalar field 
theory. 

The generalization from flat Minkowski spacetime to curved spacetime has also been 
considered. The non-relativistic Bose gas in the Einstein static universe was given by Al'taie 
10| . The extension to relativistic scalar fields was given for conformal coupling in Ref. |Tl 



and for minimal coupling in Ref. |12|. The higher dimensional version of the Einstein static 



universe was studied by Shiraishi [13]. More recently, the case of hyperbolic universes |14 



and the Taub universe have received attention. Anti-de Sitter space was studied in Ref. 
||T6| where some of the issues of our paper where considered from a different viewpoint. 

One advantage of dealing with specific spacetimes of the type mentioned above is that the 
eigenvalues of the Laplacian are known, and as a consequence the partition function and the 
thermodynamic potential can be obtained in closed form. Another approach to studying 
BEC is to try to keep the spacetime fairly general, and to calculate the thermodynamic 
potential only in the high temperature limit. This has been done by a variety of people 
||T7|^3|. In particular the symmetry breaking interpretation of BEC was given in Refs. 
2^ , |23[| . The effects of interactions have been given recently . 



The purpose of the present paper is to re-examine BEC in the case where the spatial 
manifold is compact. We will be particularly concerned with the Einstein static universe for 
which the spatial manifold is S^. Because the volume is finite the general theory presented 
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in Refs. [0,^ must be modified since the results tacitly assumed the infinite volume limit. 
Although the thermodynamics of finite volume systems has been well studied in classical 



nonrelativistic statistical mechanics p5|-p^, there appears to be some confusion in the lit- 
erature in the context of relativistic field theory, particularly in the calculation of critical 
temperatures. After a brief review of the general theory presented in Refs. P^,pB|, we will 



specialize to the Einstein static universe. The generalised ( -function |^ will be evaluated 



for this spacetime and used to obtain the effective action for a complex scalar field with 
a finite charge. Expressions for the pressure and charge are obtained. We then present a 
detailed analysis of whether BEC and symmetry breaking occur, and conclude that they 



cannot, in contrast to the expectations of Refs. p2| , p3| . The reason for the difference is 
linked to the fact that the spatial manifold has a finite volume, and we present a discussion 
of this important point. In appendix |^ we show how the generalised ( -function can be 
used to relate the thermodynamic potential directly to the effective action. In appendix 
we present a short discussion of how the analysis proceeds for the anisotropic spacetime 
obtained by identification of antipodal points on S^. 

II. THE EFFECTIVE ACTION AND BEC 

We will consider a complex scalar field defined on a (D + l)-dimensional spacetime 
manifold M ~ i? x S where S is a compact Riemannian manifold of dimension D. After a 
Wick rotation to imaginary time the configuration space action is 

r . (1,. . ,2 , 1,,. , ,,„. ^2 



fp f ( 1 1 
S[(p] = dt da^l-{ipi - W2)^ + 2 (^^2 + Wi' 



+ ^|V¥^iP + ^|V^2p + \{m' + miifl + (2.1) 

(Here |Vv?iP = g^-'di(pidj(pi with gij the Reimannian metric on S. dax is the invariant 
volume element on S.) This may be found as described by Kapusta ^ in fiat spacetime 
starting with the phase space path integral and incorporating the conserved charge using a 
Lagrange multiplier /i. The partition function may be expressed as a path integral over all 
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fields periodic in time with period (3 = T ^ where T is the temperature. 

It is straightforward to compute the effective action using the background field method 



| 35| . If we take the background field to be V5i(x) = (/^(x), ^^2(3^) = 0, then the effective action 
is 



r[(^] = 5[<^] + -lndet{/%-b]} 



(2.2) 



using condensed notation . / is a unit of length introduced to keep the argument of the 
logarithm dimensionless. 

Let {ip]^{x)} be a complete set of solutions to 



[-V^ + iR](pN{x) = (JnVn{x) 



(2.3) 



normalized by 



da^LpN{,x)LpN'{x) = 6nN' 



(2.4) 



If we define a generalized C-f unction by 



j=—oo N 



(2.5) 



then it can be shown that 23 



r[^] = %]-C^(o) + c^(o) in/2. 



(2.6) 



In order for the sum in ( |2.5| ) to be well defined, /i must be restricted by ctq + ~ /^^ > 0, 
where ctq is the smallest eigenvalue of the set. 



The background scalar field must satisfy 

6rm 



0. 



(2.7) 



The differentiation is to be computed with fi, T, V, gij all held fixed. This is equivalent to 
minimizing the Helmholtz free energy with Q, rather than /z, held fixed [^. Because C,v 
has no explicit dependence on using (|2.6| ) in (|2.7|) gives 



- V^(^ + (m^ - /i^ + ^i?) = 0. 



(2.8) 



We can expand Lp in terms of the Lp^ , which satisfy ( p.3| 



(2.9) 



where the expansion coefficients Cat are to be determined. Use of (|2.9|) in p.8|) leads to 



(Tat + m — /i )CAr = 0. 



(2.10) 



If /i^ < (To + then the only possible solution to ( p. 101) is Cat = for all A^. This leads to 
(p{x) = as the ground state and no symmetry breaking. If, however, it is possible for \i to 
reach the critical value /ic defined by 



/i^ = (Jo + 



(2.11) 



then Co in (|2.10|) is undetermined. In this case the ground state in ( p.9|) is 



(2.12) 



and there is symmetry breaking. 

In order to find the relation with BEC, consider the charge which is given by ( in units 
with e = 1) 



Q 



From (12.61) it is clear that we can write 



1 dV 



(2.13) 



Q = Qo + Qi 



(2.14) 



where 



(2.15) 



Q, = ^{cm-uo)^^i'}- 



(2.16) 



Use of the high temperature expansions in Ref. shows that for T large enough, it is 
always possible to have fi < fic and = 0. In this case Qq = 0. When the temperature 
drops, /i increases. If it is possible for /i to reach the value fic defined in ( p.ll| ), <^ 7^ and 
Qo 7^ 0. The value of T at which fi = fic is defined to be the critical temperature Tc. It was 
shown in Ref. that 



(2.17) 



exactly as in fiat spacetime (allowing for a different value of Hc)- For T < T^., we have 



Qo = Q 



-'0 



(2.18) 



T 



1 

2\ 2 



(2.19) 



These results assume that it is possible for fi to reach the critical value fic at a finite 
temperature. However without a detailed analysis it is difficult to see if this is possible. 
For example, in the non-relativistic case, when S is a fiat 2 - dimensional space, fi never 
reaches /ic and therefore there is no BEC. ( See Ref. for example.) Another case where 
fi never reaches /ic is for S a fiat 3 - dimensional space with an externally applied constant 



magnetic field More generally, if S has a finite volume extreme care must be used. In 
the past, compact spaces have been analyzed using high temperature expansions. A proper 
analysis of BEC may require temperatures outside the range where the high temperature 
expansion is valid. Use of the high temperature expansion may give misleading results as 
we show later. A recent discussion of the existence of BEC using generalised ( -functions 



has recently been given [B7 



III. THE EINSTEIN UNIVERSE 



At this point it is instructive to consider a specific example. One space which has received 
considerable attention in the past is the Einstein static universe, ie S = S''^ which has radius 
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a, with scalar curvature R = 6a~^. We will take Ui = ^-R, where ^ must satisfy the condition 
^ > |(1 — a?Tn?) but is otherwise arbitrary. High temperature expansions for the case ^ = 
where first obtained by Parker and Zhang |12[, whilest the case ^ = \ has been studied by 
Singh and Pathria [O]. Previously the critical temperature has been calculated to be 



3q 

m 



1 + 



6e 

im?a? 



(3.1) 



(see and for the case ^ = 0). q = Q /V = Q/27r^a^ is the charge density. 

It is well known [^9| -|4l|| that the eigenvalues of — on 5''^ are —N{N + 2)a~^, with a 
degeneracy of (A^+1)^, {N = 0, 1, 2, . . .). Hence = [A^(A^ + 2) + 6,^]/a^ and our generalised 
zeta function is 

(iV+l)2 



Us) 



oo oo 

E E 



(3.2) 



,=-00 to iCfj + ^^^^ + m^Y 

A slight simplification in the form of ( p.2|) is obtained by completing the square in N in 
the denominator and then relabelling the sum, n = N + 1. Thus we consider the analytic 
continuation of 



00 00 



c(^)= E E 



n 



(3.3) 



j=_^ti [{aj + ihf + a-n? + c]' ' 

where a = 2n/P, h = a = 1/ a? and c = rn? + (6^ — \)/a?. The techniques which we shall 
use are based on those of Elizalde p3|-p5| . 
We start by making use of the identity 

a-' = —— / dt f~^exp(-at), (3.4) 

r(s) Jo 

and, expanding out the (aj + ibY to get 

((s) = —— / dt f y n^exp(-anh) V exp(-aft - 2iajbt - ct + hh). (3.5) 
(We shall justify interchanging the sums with integrals shortly). Now, following the notation 



of Whittaker and Watson [42| for theta functions. 
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j2 ^-aj't-2ajbti = i^2j2 e-^' cos{2abjt) 



j=-oo 



abt 



IT 



(3.6) 



where 6*3 (2|r) is defined as {q = e'^*'^) 

00 

Osizlr) = 93{z,q) = 1 + 2 ^ g"' cos2nz. 
One immediately notes that 



(3.7) 



n=l 



00 Q 

n=l 



'(9(at) 



at \ 

TT / 



(3.8) 



An important property of the theta functions is that they are uniformly convergent for 
Xm r > 0, as is their derivative. Thus the integrand in ( p.5| ) is the product of a number 
multiphed by two uniformly convergent sums. This justifies our previous manipulations. 

We are now in a position to rewrite our zeta function in terms of theta functions, 

1 „ 1 r„ f d 



dt f-^e-^'-'"^' 



-e. 



TT 



(3.9) 



T{s) Jo ^ " \ d{at) 

Of course, although apparently much simpler, the integral as it stands could only be per- 
formed numerically since there is no representation for the theta function, other than in 
terms of the infinite sum above. However progress can be made by making use of Jacobi's 
imaginary transformation for theta functions: 



03{z\r) = {-ir) 2 exp ( — ) 03 

yiTtT J \T 

Applying this to the second theta function gives us 

1 7r5 rco 



(3.10) 



C(^) 



Tis) 



« ^0 n=l 



2-an-'t 



00 .22 

J T 



1 + 2 ^ e cosh(2j— tt) 



(3.11) 



We now split this expression into the sum of two terms as follows: 



as) = h + h, 



(3.12) 



where 



r s a 7o ^1 



(3.13a) 



and 



^, , > > cosh I 



a I Jo 



dt f 2 exp 



— (c + an )t — 



aH 



(3.13b) 



If we now differentiate ( ^.lUI ) w.r.t at , we obtain the useful identity 



oo 
n=l 



4(at)t ^ 2{at)i 



^exp 



n=l 



at 



at 



(JXJ 

exp 



n=l 



at 



(3.14) 



which, together with a standard integral representaion for the modified Bessel function 



gives us 

h = 
where 



TIC 



2-s 



ig)7-.U-exp(-.-g 



(3.15) 



y + ^ L, . E Vi^{^^-^+2(^n) - X^K^.UXn)} , (3.16) 



4(s — l)(s — 2)aa2 aQ;2r(s) „=i 



a2 



(3.17) 



Similarly, 



ar(s) 



^ cosh f ^) f; 2n2(c + an2)-^+^ 

7 = 1 \ ^ J n=l 



j-n^c + an ) 2 



where 



2vrj 



(c + ara^) 2 



i^_.+ i(^,„), (3.18) 



(3.19) 



(Note that we have not used equation ( p.l4| ) here). It can be shown that the sums in Ji and 
I2 are uniformly convergent for all s (this is because K^{z) falls off exponentially for large 
z, and = K^^). Thus we finally obtain the expansion 



TTC 



,2-s 



4(s - l)(s - 2)aQ;2 aa2T{s) 

n=l 



27r2 
aV{s] 



cosh 



/27rj6' 



£2n2(c + 



an 



n=l 



j7r(c + an ) 2 



(3.20) 



Now, for our application, we are only interested in ({0) and C'(0)- At first sight it looks 
like our effective action will be extremely complicated, however 

^ s + js^ + 0{s^), ass^O, (3.21) 



r(.) 

and so we are only left with the first term in C(0)- Similarly, when we differentiate with 
respect to s and then evaluate at s = , the result is to differentiate the first term and then 



simply to remove the gamma functions from all subsequent terms. Furthermore, using |^ 



ir.(z) = (^)%-^ (3.22) 
the double sum term in C'(0) becomes 

j=l j n=l I a j .^^-^ .^^ j 



n=l 



- ^ |ln (l - e-"^^ 

+ ln(l -e-"^"-)} (3.23b) 



where 



27r 2\ 27r6 

iUn± = — [c + an ) ± (3.23c) 

a a 

which is the usual statistical mechanical contribution. (See appendix 

Putting all of the above together, and rewriting the effective action in terms of the 
physical variables, we obtain 

r = Sm + ^Pa'c' HPc) - ^Pa'c' 

°° 1 111 

-2/3a^c^ ^ lK2{2'Knac2) - 27mac^ Ks{2irnac^)\ 

n=i {2'Rnac'i) ^ 

oo 

+ ^ nMn [(l - e-^(""+^)) (l - e-^^"""^))] (3.24) 



n=l 



where = {{n/aY + c)^ , c = + (6^ — l)/a^ and we expect S[0\ to be zero for T > Tc. 

In the above we have been exclusively concerned with the temperature zeta function. 
Originally it was hoped that by concentrating on the full sum one may be able to find 



10 



some simplification which would not be apparent if the sum was separated into its zero 
temperature and finite temperature contributions. One approach would be try to simplify 
the sum over n by using ( |3.14| ) in ( |3.13b| ). Making liberal use of the summation formulas in 
the appendix of reference one finally obtains 



where 



-c^ ln(/^c) - 



16 



32 



P i=i 3 



l=—oo 



A = 27T 



.27r/,2' 



(3.25) 



(3.26) 



y = exp(— aA) 



and the gn{y) are the generalised Bose- Einstein integrals, defined by 



9n{y) 



dx 



T{n) Jo y ^e^ — 1 

oo I 

T- 

^ In 

1=1 



(3.27) 



(3.28) 
(3.29) 



(see [0 appendix D for a discussion of the properties of these functions). Differentiating 
with respect to fi gives a conserved charge of 



^ _ 27ia^c g sinh(/?/ij) 



J 



MPc^j) - 2a' 



[c-(/. + z^)2]l 



-aA 



- 1 



(3.30) 



in agreement with |TT|], where they used the Poisson-Summation formula to derive this 
result. Examination of this result shows that it is convergent for fi^ < c; however (p.24|) is 
convergent in the larger range /i^ < c+ 1/a^. The reason for this discrepancy is most easily 
understood by consideration of the statistical mechanical term in ( |3.24| ). We can then see 
how the use of the Poisson-Summation formula results in a shortened convergence regime. 

The Poisson-Summation formula essentially replaces the original sum with a new sum 
over the Fourier modes of the original. For an even summand it becomes: 
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E/W = o E / f{n)e'-''''dn. (3.31) 

n=l ^g=-oo"'-°° 



In our case, 

27r Ox 1 ,n\ / 27r , 



/(n) = n In 



ZTT 9 1 \ / 27r r, 1 

1 - exp [(c + an ) ^ - Oj 1 ( 1 - exp [(c + an ) 2 + b\ 



(3.32) 



The analytic structure of f{n) in the complex plane is shown in Fig. |T|. There are isolated 
essential singularities at n = ±[6^ — c/a]^, and a branch cut between ±i(c/a)^. When 
< c the singularities lie on the imaginary axis and hence play no role; however for > c 
the singular points are real and tend to 1 as 6^ — c + at which point the original sum 
diverges. Use of the Poisson-Summation formula implicitly assumes that f{n) is analytic so 
that its Fourier integral decomposition exists and hence cannot be used for I? > c. Similarly, 
the approach above also implicitly assumes that the function is analytic along the real line 
in the evaluation of the Bessel function integrals. In performing the analytic continuations 
one should be careful to ensure that the domain of convergence of the final result is the same 
as that of the original zeta function. 
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FIGURES 




FIG. 1. Analytic structure of f{n) in the complex n plane, showing the position of branch cuts 
and singularities for the case > c. 



Returning now to the full result of equation (|3.24| ) we can now proceed to write down 
the thermodynamic potentials. Remembering that F = —q = —(3PV, and using the volume 
of S*^, V = 27T^a^, we can immediately write down the pressure (for T > T^) 



P= -^cHn{fc)+ ^""^ 



327r2 

^2 oo 



647r2 



C 1 1 1 1 



TT" {2TTnac 
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-| oo 

^ In [(l - e-^(""+^)) (l - g-^^'""-^))] . (3.33) 
''^ P „=i 



The first few terms in tlie pressure are essentially unobservable renormalisation constants 
(since one can only really measure differences in pressure), and the last term is the term that 
one would predict from normal statistical mechanics. The Bessel function term, however, 
is observable since it varies with volume; it is a contribution to the pressure due to the 
Casimir effect. Notice (as one would expect) that it is independent of both the temperature 
and chemical potential - its origin is in the vacuum, not the finite temperature and chemical 
potential effects. Also, because of the asymptotic form of the Bessel functions, this term 
tends towards zero as the volume (and therefore a) tends to infinity; it will not survive the 
infinite volume limit. 

Equation (|3.33|) on its own does not completely determine the pressure: /i is unknown. 



To rectify this, we must now consider the vacuum expectation value of the charge. Above 
the critical temperature it is given by ( p.l3|) 



^2 ^ 



n=l ^ n=l ^ 

(^e/3(a;„+/x) _ x)(e/3('-'n-/^) - 1) ^ 

The charge Q is made up from two parts: one contribution from particles (chemical potential 
/i), and a contribution from antiparticles (chemical potential — /i). Its high temperature 



behaviour is examined in The low temperature behaviour {(3 oo) would be trivial 
to calculate, except that a priori one can not be sure that ( |3.34| ) still holds. 

If one plots (|3.34|) as a function of /i one obtains a graph similar to Fig. |] .The charge 



diverges whenever /i tends to one of the oj^ and it is multivalued, being divided up into 
intervals by these values. Now from ( |3.34| ) /i and its' derivatives are well defined within each 
of these intervals (since the sums are absolutely convergent) hence /i is a continuous function 
of T in these regions. Furthermore, within each of these regions Q takes on all values; hence 
after choosing one of these regions one can never evolve out of it by a change of Q or T. 
Since a change of Q ^ —Q (i.e swapping particles with antiparticles) must cause /i —fi 
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only the first region is physical. 
^ 100 

50 



-50 

-100 

FIG. 2. Plot of charge Q against /x. 

Although we cannot solve for ^ analytically, in the case of the Einstein Universe we have 
done so numerically. The results are shown in Fig. |^ (which uses the charge density p = Q/V 
rather than Q. From the figure, one can see that in fact there is no critical temperature: 
fi tends towards its critical value, only reaching it at T = 0. This is most marked in the 
case a = 0.25. When the volume becomes large, the curve for fi starts to look more like 
the scenario given above and fi tends asymptotically to its critical value (but it still does 
not reach it until T = 0). This means that the S[ip] term in F is always zero : there is no 
symmetry breaking on the Einstein universe and hence BEC does not occur. 
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FIG. 3. Plots of fi against Temperature. Curves have m? = 2 and = 1/6. The dotted hne 
corresponds to a = 0.25, the dashed hne is a = 1.0 and the sohd hne is a = 20. 



IV. FINITE VOLUME SPACES 

In the previous section we examined numerically the Einstein static universe and showed 
that the critical temperature is zero, in contrast to calculations based on high temperature 
expansions (which implicitly assume that there is a non zero critical temperature!). We now 
wish to consider the generic case of a compact spatial manifold. 

Usually in statistical mechanics, we are interested in the limit as both the volume and 
the number of particles tends to infinity. It is then meaningless to talk about Q, because 
if we imagine starting off at some low temperature with no anti-particles the charge will 
also become infinite. Thus we now restrict our attention to p = Q/V. For a general finite 
volume space, this will be given by (see appendix.^ or ^) 
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_ 2 sinh(/3/i) ~ n^e^'^" 

where the ci;„ satisfy < cui < u;2 < "^3 < • • • but are otherwise arbitrary. 

The form of p vs /i must be similar to fig.l: p — oo whenever p — tu^ and by similar 

arguments to those in section |IT1| p lies in the region |p| < uji (the summand behaves like 

exp —(3uJn as a;„ becomes large - which we expect to happen as n ^ oo - therefore the sum 

is absolutely convergent). Thus we separate off the ground state contribution from p: 

2sinh(/5p)e^'^i , , , 

P = (e/3K+M)-i)(e/3K-M)-i) + (^-2) 

and we consider the situation when p is close to its critical value. We will define p = cui — 5, 
6 « 1. pe(/3, p) can be interpreted as the charge density in excited states. 

Here (|4.2| ) is to be considered as an equation in the unknown p; p is a measurable 
quantity which here we consider as fixed. Furthermore, all the equations which were used 
to derive (|4.2|) are well defined at every point except the critical value (ie are well defined 
for all 5 7^ 0); thus no matter how close to the critical point we are, we must always be able 
to solve for p. But this is just the statement that the limit of the r.h.s of ( |4.2| ) exists as 
6^0 and is equal to p. Now pe{Pc, p) tends to some limit, pei say, as 5 — > 0.(This must 
be true because the only divergent part is in the ground state; the sum in pe is absolutely 
convergent and therefore must sum to some value at p = cji). Let us call the ground state 
contribution in ( [4.2| ) pg(/3,p) and let us also define the constant pi = p — pei. 

Assume that we have a critical point, ie a point where p — * cji as /5 — (3c- Thus we must 
have that 

Pg{P,uj^-5)^ pi as5-.0. (4.3) 
But by expanding pg in terms of 5, it is easy to show that 

^ o(l-\ as 5-^0 (4.4) 



(35V W, 

The only way that we can get a finite limit for p as 5 ^ (and (3 (3c) is if (3V = 0(1/5) 
as 5 — > 0. 
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For a finite volume space we must have P = 0{1/S) as S ^ 0; thus P ^ as fi ^ ui 
and the critical temperature is zero. Since for all finite (3, fi ^ uji this also means that there 
can be no symmetry breaking in such spaces. However, as we saw in the Einstein universe, 
if the volume is large we expect that fi can be quite close to its critical value for some finite 
temperature; thus there will be a large number of particles in the ground state and the 
physics will look similar to the infinite volume physics. In this situation the temperature at 
which we expect some degree of condensation is at a critical temperature obtained by taking 
the infinite volume limit. Alternatively one can try to define the onset of condensation by 
looking for the maximum of the specific heat since this is the property which is usually 
searched for experimentally. 



We examined the Einstein static universe numerically, and showed that the critical tem- 
perature defined by /i — > tui is zero in contrary to previous calculations. These gave the 
wrong result because they were based on high temperature expansions. We then generalized 
our result to arbitrary curved spaces of finite volume and showed that there is no symmetry 
breaking on these spaces at finite temperature and density. 
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APPENDIX A: NON INTERACTING SCALAR FIELD AT FINITE 

TEMPERATURE 

In this appendix we wish to show in detail how the generalised zeta function techniques 
reproduce the earlier results of 0. Taking the field theory of section. |I| we define {(pn} s.t. 



V. CONCLUSION 



ACKNOWLEDGMENTS 




(Al) 



18 



If we let a = 27t/P then 



(A2) 



n j=-co 

where the sum over n in ( [A2| ) is over all the eigenvalues of ( |A1| ) . Using ( |3.4| ) and the Jacobi 
identity we obtain 



t-l (Eexpl-.?,)) (l + 2|:exp cosh (^)\ (A3) 



aTis 



The integral in (|A3|) naturally splits up into two parts. In the second part, the integral is 
easily done using the Bessel function integral (by a change of variable t — > elt): 



2712 

aV{s) 



OO n / 

EE/ rftf-iexp 

n 1=1 "'0 \ 



aH 



cosh 



(^) 



47r2 



E E cosh 



f2TTjfl\ fTCj 



\ a J \ a 



'2 1. 



/27rjen 
\ a 



Differentiating this, then evaluating at s = we obtain (remembering ( 3.21|) 



(A4) 



a 



n j=l 



\ a J \TTj 



2 V a 



n j=l J 

-E{Mi 



,-/3(en+M)' 



+ ln(l-e-^(""+'^))} (A5) 



where we have used the fact that /? = 2n/a and ln(l — x) = — J2'^=i x""- 
The first integral in ( |A^) is 



7r2 



ar(s) „ Jo 



E / ^ exp(-£^t) = ^ ^ 



a r(s) 

7ri r(s-|) 
a T(s) 



p2s-l 
n 



Q(s-l/2) 



(A6) 



We expect Ce^ (s — 1/2) to be an analytic function of s except for a few simple poles. This 
implies that it has a well defined Laurent expansion about s = 0, which we will write as 



(el-^ + Ao + 2A,s + 0{s'). 



(AT) 
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Note that A^i is the residue of the generahsed zeta function, using as the eigenvalues e^, at 
the point s = 1/2. This expansion, together with ( |3.21D , immediately gives 



7r2 



T(s) 



s=0 



Ti^T I -- ) A_i = -2nA_i. 



(A8) 



Also 



TC2 — L 

Tis) 



CeAs-l/2) 



+ 7r2r(s- -)(Ao + 7) + ass^O, 



(A9) 



where ip{z) = r'(z)/r(z). Making use of ■?/'(— 1/2) = —7 + 2(1 — In 2) and then substituting 
in (|2.6|) (remembering that the S[i^] term never appears - see section |IV|) gives 



(AlO) 



B, we 



r = pA.i In /2 + /3{(2 - 2 In 2 - 7) + + 7} 

+ In [(1 - e-^(^"+^)) (1 - e-^(^"-'') 

By comparison with the result of Haber and Weldon P|, or from examination of 
note that the e are the energy eigenvalues and identify 

^ = A_i In /2 + { (2 - 2 In 2 - 7)^.1 + + 7} 

n 

which is essentially the (regularized) Casimir energy. 
Now , because = cr„ + m^, 

1 



(All) 



dt t^-ie-*('^"+™ ) 



T(s) JO 

It may be noted that G(t) is tr exp[— 1(— + Ui)] which has a known asymptotic expansion 



-f- / dt f-^e-™ *e(t). 

(s) Jo 



(A12) 



as if: — s> namely. 



e(t) ^ (47rt)-^/2 Yl t'(^k. 

k=0,l/2,l,... 



(A13) 



The 9k are the same coefficients that are used in , and depend only on the geometry and 
the conditions on the fields at the boundaries of S. (Typically they involve products of the 
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curvature invariants and the masses of the fields on S.) D is the dimension of S and this 

is the only place in which it enters explicitly. Use of this expansion in ( |A12| ) together with 

expanding the exp(— m^t) allows one to find 

47r-^/2^£{?(-l)^/2-'ei/2+im^-2' for D even 
47r-^/2 j:l^;l'^^/\-iyD+i)/2-i0^ m^+l-2l for D odd ^^^^^ 

Thus the residue of Ce,i(0) exists and can be written down in general. As claimed earlier, 
one can see that it depends only on the geometry of the space. Unfortunately the authors 
are aware of no similar method for calculating Aq in general, and one is reduced to working 
out the analytic continuation of Ce„ case by case. However, since the statistical mechanical 
contribution can be written down explicitly, one only needs to focus on the zero temperature 
zeta functions. 



APPENDIX B: ANTIPODAL IDENTIFICATION 

A simple modification of the Einstein static universe is to identify antipodal points on 
S"^. Two cases are possible: periodic identification and antiperiodic identification of the 
fields. In the case of periodic identification, all the odd modes of — must vanish. Thus 
the generalised zeta function in this case is 

{2N + 2y 







3,N 



oo oo 

E Et 



-oo 7V=0 



(2Af+l)(2Af+3) I 6^ 



+ ^ + m' 



(Bl) 



At first sight, it appears that we shall have to repeat the entire calculation for the zeta 
function above. However if one makes the change of variable, n = + 1, and notes that 
(2A^ + 1)(2A^ + 3) = {2n - l)(2n + 1) = An^ - 1, we can rewrite (0) as 

oo oo 







E E 



4x E E 



n 



(B2) 



j=-oon=l + + + C] 

where a = 27r//3, a = A/a'^, b = fi and c = (6^ — l)/a^. This is the zeta function that we 
had before, and 
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r, = l/?aVln(/2c)--/?aV 

°° 1 111 

—f3a^c'^ ^2 \K2{7mac^) — nnac^ K3{TTnac2)\ 

n=i (irnac^) 

oo 

+4 J2 In [(l - e-^('^"+^)) (l - e-^('^"-'^))] (B3) 

n=l 

ojn = {{2n/aY + c) and c is as before. The factor of four in C, appears simply because the 
space has half the volume that it had originally .As the volume becomes large we expect the 
log term in F to also be directly proportional to the volume, which explains the four here 
(as the sum tends to an integral, the measure will give a term proportional to a^. But in 
( P3|) a a/2 compared to our original effective action, leaving us with an overall factor of 
two). The spacetime gives us an effective thermodynamic mass, c, identical to the value on 
S*^, however the detailed dependence of the sums at small volume is now radically different. 

Because the expression ( p.3| ) is an absolutely convergent double sum for all points at 
which it is defined, we can write 



Codd — Cs'i — Ceven (B4) 

which will be true even after analytic continuation. 

Previously, it has been remarked that there would be many interesting properties of this 
space associated with the possibility of have a x dependent charge density. However, this 
charge density appeared through the S[0\ term which we have seen will be absent. 
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